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ABSTRACT 


The objectives of the project described in this 
thesis were to ; 1) provide the documentation that is 
needed for a Naval Postgraduate School student to use 
the general purpose finite element computer progran 
called SAP IV, and 2) to make available, and prepare 
the users manual for, a pre-and postprocessor progran 
called SUEROUTINE PSAP. This subroutine, which was 
developed at the NASA Langley Research Center, has 
been mnodified to specifically plot the finite element 
model gecmetry for SAP IV models and to postprocess 
displacement data for those models on the NPS Calcomp 
Model 765 Plotter. The input and output for SAP IV and 

SUBROUTINE PSAP are discussed in detail. The codes 
have been used successfully in AE 4102, Flight Vehicle 
Struetural Analysis IT. 





TABLE OF CONTENTS 


308 INTRODUCTION eevsrteseeseseeseseses8te82828 @eeesese@*see@ensvsesenseeseespeeseseseeeses 8 
A. STRUCTURAL ANALYSIS SOFTWARE PACKAGES........ 8 
B. PRE-AND POSTPROCESSORS FOR STRUCTURAL ANALYSIS 


PROGRAMS .- eee ee ee a ee ee ee ee cece sces Sree 5 oe 3 
Ce theSiS MOTIVATION. ....<-- erence ccces cece ces es 2 
erect Re tOebHE USE GF SAP LTV AT NPS. cescewecnacess 11 
A. DESCRIPTION OP SAP [Vewceeeneveeee re eee 

B. COMPUTER CARD DECK PREPARATION AT NPS FOR 


SAP RIAU PewemeweteWat eter che ts «oa: cele el cicokere! eel slcisleels ee 6 ela es 66 6 64 6 #060 6% 12 


Ge oELPPUL POTNTS ON DATA DECK PREPARATION. .«s..- 15 

Dee Geert LONMOr SAP IV RT NPS sce se ee cee ees we 16 

iif. GUIDE TO THE USE OF SUBROUTINE PSAP, A PRE-AND 
ONO Cho reenOR GAL (vials ne ais se ees kicecaannssweacecsenas 20 
De SUSROUL PN E VUES CREPTION Sec cca ce cs wees wna wae ee «820 
Reese oOUleN Perot er Uo Bela sle cals sewececee ses aecisa 21 
jeeeecenecal Sct-wp Of Input DecK.sccsacansses 21 
PamwOOPCORCEOU CabdS. 6.6. cence cece wanes necens «622 

Seo blebs TEOCISCiNG Gan 5 obo po Gon ene eee 65) 


Ure See ara do to never aie, aan 6 eee) ceca Bye ers a ele Si ciedecasats so SEO 
PRN amer DS Ce OPT WON ecco cise «eee iss Gs es ere Ale) 
So (GreumSicies, Qcoyewles Waves! (Eek 4 oye eer ool 28 
fempeasec INOeiriracat1on, Card... 2. <2... 26 desea sal 
8. Deck of Displacement Data to be Plotted.. 33 
So WB REIL use EI UCTE S SS Gea no Meee a Geers e esac 4 


See ET HOD POR VALTER ENG SUBROUTINE PSAP.......... 38 
Dee StGNEPECANT ASPECTS OF SUBROUTINE PSAP....... 39 
leeeexoleded Plots... 5. SoGnbe SOR nen ee me are Ta wo) 
Ce OmmnOno mon NOGCES sense sese saivecccewecesec 19 
Swope GiLlcatiOn GCL VieWeu cass esecccerscesss 40 


POSSI DUCmSOURGES Ome Libt Ol Sic oes slele scl s'si ele) os 40 





ivemreceanG LUSTONS AND RECOMMENDATIONS. ....cnceccesceses 42 
mependix At EXAMPLE PROBLEM USING SAP IV. .ssessacscceee 44 
Appendix Bz: EXAMPLE PROBLEM USING SUBROUTINE PSAD...... 46 
Appendix C: SAMPLE SUBROUTINE PSAP OUTPUTS FOR A FINITE 
BEEMENT MODEL OF A WING. cc ce cccncccccccscccsccsccsecsess § 33 


Mependix D: SUBROUTINE PSAP LISTING........ Sere ayes ee 209 
IG t Se OF FIGURES. ees eee eo scee#ses#2e2e8 oe eeeeses#ssseeseesee#et82f 8288 eeseesee@ess#ses#ss#se8e¢e 60 
JE IOS BE OF REFERENCES. eeseesedes# @ @eeoeees#*e?d@eeseee@eesee@ess#ensteet eee @®eeee#s#3es¢6 89 


INITIAL DISTRIBUTION ONES Ul emeMeMeMetel site el oreleleleccterelel ote a sie «sec icrecere 90 





ACKNOWLEDGEMENT 


I wish to acknowledge the help that I have received in 
the course of my research. 

I am particularly indebted to Professor R. B. Ball, who 
provided me with the guidance and aid necessary to complete 
the project. 

I am also grateful to Professor Gilles Cantin, who 
provided me with valuable information in starting and 
carrying through my work. 

I also wish to give credit to the NPS programming 
consultants, and in particular Sharon D. Raney, for the fine 
assistance I received in debugging my computer runs. 

My gratitude is also heartfully extended to Richard 
Citerley for providing me with the computer card decks 
necessary to begin my project. 

I especially want to express my Sincere appreciation to 
my wife, Karen, for her patience and perseverance with the 


typing of this thesis. 





I. INTRODUCTION 


See ee ee 


A. STRUCTURAL ANALYSIS SOFTWARE PACKAGES 


The past decade has seen great strides in the field of 
computer software packages that have been developed fcr use 
in structural analysis. A great deal of time and effort has 
been expended in the development and distribution of these 
packages. Today there are structural analysis computer 
programs for almost every conceivable structure an engineer 
could desire to analyze. The types of analyses performed by 
a given package vary widely and may include features for 
linear or nonlinear materials, static analysis, dynamic 
analysis, buckling analysis, and nonlinear dynamic analysis, 
to name only a few. These programs can be grouped into two 


major categories - special purpose or general purpose. 


In order to expose the students at the Naval 
Postgraduate School to the use of a general purpose 
Structural analysis program, as well as to provide the 
capability of using such a program in research work, the 
finite element structural analysis program SAP IV (Reference 
1) was acquired and made operational at the Naval 
Postgraduate School by Professor Gilles Cantin of the 
Mechanical Engineering Department. SAP IV can perforn 
linear static and dynamic analyses on one-, two-, and 


three-dimensional structures. 





B. PRE-AND POSTPROCESSORS FOR STRUCTURAL ANALYSIS PROGRAMS 


After the development of a large number of structural 
analysis programs, users began to recognize that a disparity 
existed between efficient general purpose structural 
analysis programs and optimum utilization of these programs. 
Many of today's software structural analysis packages, such 
as SAP IV, require the user to prepare and reduce tremendous 
amounts of data. The need existed for some aids in 
processing and reducing these large quantities of data. 
Conseguently, there have been many pre- and postprocessors 
developed for a specific use as well as for a general use 
basis during the past several years. The value of a given 
processing package lies in its ability to aid the user in 
preparing his model geometry, in his data checks, and in 
processing the output in an easily understood fashion. One 
of the most effective means utilized in preparing or 
reducing data is through the use of visual displays, whether 
they be designed primarily for graphic presentations, such 
as the Strcmberg-Carlson, or for paper plots, such as 
Calcomp. This ability to visually display input and output 
data is a highly valuable tool for the structural analyst. 


Ce THESIS MOTIVATION 


The desire to have a pre- and postprocessor that could 
be used in conjunction with the general purpose structural 
analysis program SAP IVY prompted the acquisition and 
implementaticn of a general use plotting package by this 
author. After researching the possible options, an existing 


program was obtained from Anamet Laboratories, San Carlos, 





California. The program was originally developed at the 
Langley Research Center, Hampton, Virginia by Gary L. Giles 
for use with modern digital computers. The program, details 
of which can be found in Reference 2, generates oblique 
orthographic projections of three-dimensional finite element 
models and is distinguished by its provisions for 
generality, ease of use, different display options, and 
computational speed. The computer code was written for use 
on CDC 6000 series machines and had to be modified somewhat 
for use on the NPS 360/67. The modified version of the 
program is now available for use with the NPS Calcomp 
Plotter Mcdel 765. 


The primary purpose of this thesis is to provide the 
necessary documentation in order that students enrolled at 
the Naval Pestgraduate School, and specifically in the 
courses AE 4101, 4102 (Flight Vehicle Structural Analysis I 
and II respectively ), can, with a minimum of difficulty, 
effectively utilize SAP IV and its now-functional pre-and 
postprocessor PSAP. The remainder of this thesis is broken 
down into two major subdivisions, 

LE ESUIDE SONTH® US OF SAP Iv AT THE NAVAL 
POSTGRADUATE SCHOOL 

2 elite GULDE TO Ta 
PRE-AND POSTPROCESSOR FOR SA 
Appendices A and B of this work provide detailed examples 


B OF SUBROUNTINE PSAP, 


for the input preparation and output reduction of data using 
both SAP IV and SUBROUTINE PSAP. 
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Fi.. GUIDE TO TRE US 


A. DESCRIPTION OF SAP IV 


SAP IV is a general purpose structural analysis digital 
computer program that can provide a finite element solution 
for both the static and dynamic analysis of linear 
structural systems. aA detailed user's manual is contained 
in Reference 1. The program has the capacity to analyze very 
large three-dimensional systems, as well as small systems, 
Meta no lass in efficiency. SAP IV, which is coded in 
FORTRAN IV , is a very flexible program and can be 
considered a very efficient aid to the analyst. The purpose 
of this section is to provide the necessary additional 
documentation, above that provided in Reference 1, for a 
Student at the Naval Postgraduate School to make use of the 


program. 


The methods of analysis and the construction of the 
program are not included in this section, but can be found 
in Reference 1. The program contains nine finite elements of 
the following types: 

{a) three-dimensional truss element, 

(b) three-dimensional beam element, 

(c) plane stress and plane strain element, 

(d) two-dimensional solid element, 

(e) three- dimensional solid element, 

(f) variable-number-nodes thick shell and 
three-dimensional element, 


(g) thin plate or thin shell element, 
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(hb) boundary element, 

(1) pipe element (tangent and bend). 
There are numerous options and combinations of static and 
dynamic analysis that are available to the user of the 
program. Reference 1 provides specific details of the many 


available user options. 


B. COMPUTER CARD DECK PREPARATION AT NPS FOR SAP IV 


Figure 1 outlines the overall computer card deck 
necessary to access and utilize SAP IV as it is currently 
operational at NPS. A complete detailed breakdown of the 


necessary IBM job control cards follows on the next page. 
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4// { STANDARD GREEN JO3 CARD ) 

4//GO EXEC PGM=SAP,REGION=260K 

//STEPLIB DOD UNIT=2321,VOL=SER=CELO02,DISP=SHR, 
// DSN=FO559.SAPLM 

//GO.FTCLFOOL DD UNIT=SYSDA,SPACE=(CYL1(616)), 
// OC B=(RECFM=VS8,LRECL=1284,BLKSIZE=2572) 
//GO«FTOZFO01 OD UNIT=SYSDA,SPACE=(CYL2(675)), 
// CC B=(RECFM=VSB,LRECL=1284,BLKSIZE=2572) 
//GO-FTO3FO0OL, DD UNIT=SYSDA,SPACE=(CYL,(676)), 
// OC B= (REC FM=VSB, LRECL=1284, BLKSIZE=2572) 
//GO.FTC4FOOL DD UNIT=SYSDAySPACE=(CYL,(6,6));, 
// DC B=(RECFM=VSB,LRECL=1284, BLKSIZE=2572) 
//GO«FTCS5FOOL DD ODNAME=SYSIN 


//GO.FTC6FOO!L DD SYSOUTH=A,SPACE=(CYL, (351)),5 
el OC B=(RECFM=FBA,LRECL=133 ,BLKSIZE=1330) 
//GOeFTCTFOOL DD UNIT=SYSDA,SPACE=(CYL1(556)); 
// OC B= (REC FM=VS8, LRECL=1284, BLKSIZE=2572) 
//GO-FTC8FOOL DD UNIT=SYSDA,SPACE=(CYL1(676)), 
// OC B= (RECFM=VS8,LRECL=1284, BLKSIZE=2572) 
//GOsFTCSFOOL OD UNIT=SYSD4A,SPACE=(CYL1(6,6)), 
vf DC B=(RECFM=VSB,LRECL=1284, BLKSIZE=2572) 
//GO«FTLOFOOL DD UNIT=SYSDA,SPACE=(CYL (676) )y 
7/ OC B=(RECFM=VS8,LRECLHL284, BLKSIZE=2572) 
//GQ-FTLIFOOL DD SYSOUT=BySPACE=( TSK, (2052))4 
// DC B=(RECFM=FB,LRECL=80, 8LKSI ZE=7200) 
PAGOsoYSIN DD * 


( SAP IV --- DATA ) 


/ ( STANDARC NPS EOF CARD } 
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After the //GO.SYSIN DD * card, the deck of SAP IV data that 
is prepared according to the instructions in the appendices 
of Reference 1 follows. Following the SAP IV data is the 
standard NPS end-of-file card (/*). 


Cc. HELPFUL POINTS ON DATA DECK PREPARATION 


There are several possible areas where a user can err in 
his preparation of data for SAP IV. The following items are 
representative of a few common sources of error. 

1) Particular attention should be directed toward 
using the correct formats and correct card columns in key 
punching data on cards. (Uece, wmteqer formats, Trignt 
justiried) 

2) The program has several internal data generation 
features inherent in it. Simply because data were generated 
during the program execution does not necessarily mean they 
were generated correctly. Errors in the user- prepared input 
cards used in data generation can cause severe discrepancies 
to occur during program execution. The generated data 
should be carefully checked for accuracy. 

3) In order to terminate a given problem a number of 
blank cards are necessary at the end of the input data deck. 
Specific points of interest pertaining to the problem of 
termination can be found in Reference 1 under note (1) on 


page II.1 and in note fa on the top of page V.2. 


SAP IV has a wide range of options where a clever user 
can fully exploit the full capabilities of the program. The 
detailed decription of these features is found in Reference 
1; however, the following list summarizes a few of the nore 


useful options and features. 
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1) Data check only mode of execution. 
2) Nodal point and element data generation 
capabilities. 
3) Five different types of analysis. 
a. static analysis 
b. eigenvalue /vector solution 
c. forced dynamic response by mode superposition 
d. response spectrum analysis 
e. direct step by step integration 
4) Automatic punched computer card output of 
displacement results. This feature was added to the progran 


by this authcr to make possible graphic postprocessing. 


The foregoing discussion is not intended to be a 
complete diagnostic summary of SAP IV, but rather an aid to 
the student who desires to get started using the program as 
it currently exists on the NPS IBM 360/67. The above 
discussion, along with a copy of Reference 1, should provide 
a jumping off place for a novice structural analyst. A 
complete example of the job control cards and input data 


deck for a static truss analysis can be found in Appendix A. 


D. ALTERATICN OF SAP IV AT NPS 


PectiscmmehOuNEi is 4 vyery flexible program, it is 
possible for a user to make modifications to the basic 
program. The computer code is complex, but it is not overly 
difficult to modify parts of the program in order to satisfy 
a specific user need. A method that has proven successful 
for this author, in the nodification of a SAP IV subroutine 
to provide punched output of displacement data, and to 
create a personal version of SAP IV, is outlined in this 


section. 
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The first step is to define the subroutine name that the 
user wishes to modify. A listing of the source program is 
available through Professor R. E. Ball, Department of 
Aeronautics, OE Professor G. Cantin, Department of 
Mechanical Engineering. Having defined those portions of the 
program for which modification is desired, the next step is 
to obtain a punched copy of those desired routines by using 
the following format. 

// (STANDARD NPS JOB CARD) 

//SYSPRINT DD SYSOUT=A,SPACE=(TRK, (10,1)) 

//STSUT1 DD DISP=SHR, UNIT=2321, VOL=SER=CELOO1, 

// DSN=F0099.SAPSR 

//SYSUT2 DD SYSOUT=B 

//SYSIN DD * 
PUNCH TYPORBG=P0, MAXNAME=2 
MEMBER NAME=(SUBROUTINE NAME) 

/* 

NOTE: SUBROUTINE NAME is the name of the SAP IV subroutine 


reguired. 


With the desired routine decks now in hand, the 
necessary changes can be incorporated into the subroutine 
deck and an cbject deck is then obtained as follows: 

// (STANDARD NPS JOB CARD) 

// EXEC FORTCD 

// FORT.SYSIN DD * 

(HODIFIED FORTRAN SOURCE DECK) 


/* 


The next and final step is to take the object deck 
obtained in the previous step and insert it into the proper 
position in the control cards that are illustrated on the 


following page. 
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There are many options available in the creation of load 
module libraries and they are discussed in detail in 
Reference 4, sections II. and {iIf. The portion of the 
preceding example control card deck that would necessitate 
modification is the //ULINK.SYSLMOD Card, where 
FO559.SAPLM (SAP) should be changed according the following 
format : 

General format: Lnnnn. anyname 
where Lis : 
S-student data set 
F-faculty data set 
Qnonn is : 
user number assigned 
anyname is : 
any unique name assigned by the user (1-6 
characters in length with the first character alphabetic), 
and CELO02, should reflect an appropriate data cell with 
available space. 


The procedure outlined above will create a modified 
version of SAP IV on a chosen data cell available to the 
user for a period of 90 days. The program can now he 
executed with appropriate modifications to the //STEPLIB 
card as discussed previously in section II. B. Any questions 
concerning the creation of load module libraries or their 
execution can be answered by any of the programing 
consultants on the first floor of Ingersoll Hall. 
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A. SUBROUTINE DESCRIPTION 


SUBROUTINE PSAP is a modified version of the oblique 
orthographic projection program that is found in Appendix B 
of Reference 2. The program, originally developed for use 
at the NASA Langley Research Center, Hampton, Virginia, 
required some changes so that it could be used in 
conjunction with SAP IV and the NPS Model 765 Calcomp 
Plotter. The original version of the plotting package 
allowed for various geometry and displacement data input 
options. However, the subroutine, as currently filed, is 
constructed strictly for use with the input and output of 
SAP IV. With some slight aodifications to the subroutine, it 
could be adapted to any number of different types of input 
geometry or displacement data decks. A method EOL 


aeeompliishing this is outlined in section III. Cc. 


The current capability for generating oblique 
orthographic ‘projections of SAP IV finite element models is 
limited to the following types of elements: 

1) Type 1, Three-Dimensional Truss Elements 

2) Type 2, Three-Dimensional Bean Elements 

3) Type 3, Plane Stress Membrane Elements 

4) Type 4, Two-Dimensional Finite Elements 

5) Type 6, Plate and Shell Elements (Quadrilateral). 
The undeformed topology of the finite element model, useful 
in checking input data, as well as the displaced topology 
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projection of the same model, can be obtained from PSAP. The 
subroutine contains many different options and permutations 
of those cptions, some of which are listed below. 
1) plots of models annotated with grid point numbers 
of element numbers 
2) plots of portions of nodels 
3) exploded plots of model sections (i.e., line 
elements coincident with the edges of triangular or 
quadrilateral elements may be difficult to single 
out, program provides a capability to separate 
elements so that their absence or presence is easily 
detectable) 
4) displacements superimposed on grid point 
coordinates of the undeformed structure 
5) displacements represented as vectors extending the 
undisplaced grid points. 
There are many more combinations of options that are 
available and will be more specifically outlined in the 
remainder of this section. Appendix B details a complete 


input-output example of a pin-jointed truss. 


Be SUBROUTINE PSAP USE 


1. General Set-up of Input Deck 


In general, the correct sequence of computer cards 
required to utilize SUBROUTINE PSAP is shown schematically 
in Figure 2 and consists of eight separate major groups as 
follows: 

1) a group of JCL cards and main program to allocate 
storage 
2) a single card containing title information 


3) Namelist OPTION containing values to determine if 
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proper storage allocation is available and specifying 
various program options 
4) a geometry deck (SAP IV data deck) containing grid 
points and connectivity of the model 
5) an optional single title card used to identify the 
deck of displacement data to be plotted 
6) a single card containing the value of the total 
number of SAP IV load cases and an optional scale 
factor 
7) the deck of displacement data to be plotted 
(output of the execution of SAP IV-static analysis) 
8) Namelist PICT containing values to specify the 
type of plot desired and what information is to be 
included on the plots. 
BY repeating parts of the basic input data to the prograan, 
different plots of the same data can be generated. 





The sequence of cards on the following page depicts 
the necessary JCL to execute SUBROUTINE PSAP. 


22 





// (STANCARD GREEN JOB CARD) 
// &XEC FORTCLGP,REGION.GO=180K 


PeogeRissYoiIN CO * 


( MAIN PROGRAM ) 


/ ( STANCAREC NPS EOF CARD ) 


//LINKeUSOCD CO UNIT=2321,VOL=SER=CELOOZ,D0ISP=SHR, 


Py SUSN=F0559.PSAPLM 


//LINKeSYSIN CD * 


INCLUDE USDD( PSAP) 


ENTRY MAIN 


/ ( STANCARC NPS EQF CARD ) 


//GOoFT1OFOOL OD UNIT=SYSDA, 


J/ SPACE=(CYL4(351))5 


// SCB=(RECFM=VS,BLKSiZE=3520) 


//GOeSYSIN DO * 


eA epee oak niNG WITH THe TITLE CARD --~ S&c FIGURE 2) 


7 ( STANCARC NPS EQF CARD ) 
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/* 
REMAINDER OF INPUT DECK 
REPEATS PRRTS DEPENDING 
CN THE VALUE OF KOve iN 
NAMELIST PICT 






(FROM SAP IY 
EXECUTION? 






CARD 


7 DISPLACENENT DATA DECK 


OMIT IF 


ances? NILC, SCALEF 


GHEE 1D CRRD 


GEOMETRY INPUT DATA DECK 


NAMELIST OPTION 


(ror SAP IY) 


uCL CARDS AND MAIN 
CALLING PROGRAM 





Figure 2 > PSAP INPUT CARD SEQUENCE 
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di. ~Sarn Program 


The main program consists of three cards. ae 
allocates the proper storage and calls SUBROUTINE PSAP. An 
example of a main program for a finite element model that 
contains 400 nodes is : 

DIMENSION Z2Z (2800), DISPD (5,3,400) 
CALL PSAP (Z2ZZ, 2800,DISPD,400) 
END 
The definition of the arguments used in calling SUBROUTINE 
PSAP (222, NZ, DISPD, NON) are : 
ZZZ-blankcommon array used to store nodal coordinates 
and displacements 
NZ-lenogth of array Z2Z2Z (NZ is determined by multiplying 
the number of nodes in the model by seven, i.e., NZ=no. of 
nodes * 7) 
DISPD-a three dimensional working array used in 
subroutine DATA9 : 
NON-number of nodes in the model 
It is crucial to dimension Z22ZZ (7*NON) and DISPD (5,3,NON) 
correctly in the main progran. Improper dimensioning of 
these arrays can cause output errors that are not readily 


traceable. 


Gee litle Card 


This Card contains any desired alphanumeric 
information in columns 1 to 80. The title will appear at the 
beginning of the plots. 
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5. Namelist OPTION 


=< SSeS Se 


Namelists are a convenient means of inputing the 
names of several parameters along with their corresponding 
values. For the NPS IBM 360/67, the format for using 
namelists is as follows: 

1) card 1-&name-beginning in card column 2 where name 

is the name of the subject nanelist 

2) succeeding cards- beginning in card column 2, the 

Mames of the variables and their values separated by 
4 commas 

3) final card -&END-starting in card column 2. 


An example of a namelist format is shown below: 


col 

&OPTION 

NNDEST = 400,NUDISP =1, 

PSIZE =8.0 

&END 
Any or all parts of a defined namelist may be included, and 
each parameter may be specified in any order between the 
6name card and the &END card. 


The description of the variable names in Namelist 
OPTION and their default values are contained in Reference 
2. They are given here to assist the user in data 


preparation. 


[aa eo a [SS ee ae ee See 


NNDEST - 1 The number of nodes(NON) as defined in tae 
program 
NBDISP - 1 0 x~direction displacements not to be input 


1 x-direction displacements to be input 
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NOLS? = 


NWDISP - 


KGEOM - 


mala — 9 


9 


1 


1 


NVALUS - 0G 


meee SRO — 


KELOCT - 


BorackE = 


1 


1 


316 10) 


y-direction displacements not to be input 
y-direction displacements to be input 


z-direction displacements not to be input 


a © mm © 


z-direction displacements to be input 

( NOTE: when SAP IV displacement data is to 
be used, NUDISP=NVDISP=NWDISP=1; for no 
displacement data NUDISP=NVDISP=NWDISP=0) 
Specifies the subroutine and corresponding 
method of input for model geometry 

1 subroutine GEOM1, a user-supplied 
subroutine 

2 subroutine GEOM2, a user-supplied 
subroutine 

9 subroutine GEOM9, reads in grid points and 
element data specifically from a SAP IV data 
deck 

Specifies the subroutine and corresponding 
method of input for displacement data 

1 subroutine DATA1, a user-supplied 
subroutine 

2 subroutine DATA2, a user-supplied 
subroutine 

9 subroutine DATA9, reads a punched output 
displacement deck from execution of SAP IV 
NOT INCORPORATED-ALLOW DEFAULT 

Grid point numbers are stored in the proegraa 
from 1 to the total number of grid points 

Q no internal resequencing of grid points 
necessary; they are already in ascending 
order starting with 1 

1 resequence grid points from lowest grid 
point number to highest grid point number 
Specifies the type of output device to be 
used ( ALLOW DEFAULT ) 

Space between plots in the y-direction,in 


inches( see Figure 3 for an explanation of 
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axis orientation ) 

Potse = 10.0 Paper size in x-direction, in inches( used in 
scaling of plots to insure this dimension is 
not exceeded) 

MOCASE ~- 0 Q no identification card preceeds the deck 
of displacement values 
1 identification card preceeds the deck 


of displacement values 


6. Geometry Input Data Deck 


—2 =o 


This portion of the input deck contains the grid 
point locations and the element connectivity. The deck has 
one of the following forms, depending on the value of KGEOM 
in the Namelist OPTION. 

KGEQM = 9 
AGEUR = 3 

Calls subroutine GEOM9, which is constructed to read 
SAP IV geometry data. 


(a) When KGEOM is specified as 9, the complete input 
deck of SAP IV data, prepared according to Appendices I 
through IV of Reference 1, is placed after the 6END card of 
Namelist OPTION. The portions of SAP IV data deck that 
involve load, mass, or dynamic analysis data are not part of 
the input geometry data to SUBROUTINE PSAP. Only the grid 
point locations and the element connectivity data are used 


to generate the orthographic projection of the model. 
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(b) It is possible to input only a portion of the 
finite element model for a data check. To do so, it is 
necessary te modify SAP IV element control cards (described 
in detail in Reference 1, Appendix IV) for the desired 
element types to reflect the portion of the element 
connectivity cards that will be input. For example, if only 
connectivity for element numbers 15 through 50 of element 
type 1 (truss elements) of a SAP IV data deck are available 
for input, it is necessary to alter the element control card 
for the truss elements. Ihe field on the SAP IV control card 
that defines the total number of truss elements (card 
columns 6-10) would reflect the upper bound, in this 
case-50, and card columns 65-70 would reflect the lower 
bound, in this case-15. All nodal coordinates for the 
entire mcdel may be input, or only those that specifically 
define the portion of the finite element model to he 
plotted. In either case, the nodal coordinates that relate 
to element numbers 15-50 must be specified. Unknown results 
Will occur when trying to plot an element whose node points 
are not specified. The above feature is valuable in a case 
where several different people are preparing different parts 
of a large data base for a SAP IV problem and desire to 
individually check their inputs grapnically for accuracy. 
KGEOM = 1 

Calls subroutine GEOM1, which is prepared by the user 
to read geometry data from a program other than SAP IV. 
KGEOM = 2 

Calls subroutine GEOM2, which iS prepared by the user 


to read geometry data from a program other than SAP IV. 


Use of KGEOM=1 or 2 requires modification of 
SUBROUTINE PSAP to fit the specific format of the user's 
input gecmetry data. A method for doing this will be 
discussed in paragarph C. of this section. 
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ndentification Card 


If IDCASE =0 is specified in Namelist OPTION, this 
card is omitted. The card, if present, contains any desired 
alphanumeric information in card columns 1-80 which will 
identify all displacements for a given case. For IDCASE = 1 
and SAP IV punched displacement data, a case identification 
card must appear before each Namelist PICL for every 
different load case that is plotted in addition to load case 
one. This is illustrated in Figure 4. A maximum number of 
five different load cases can be obtained from SAP IV. The 
case ID card information, if present, will appear before 
each load case's DISPLACEMENT DATA TO BE PLOTTED section in 
the printed computer output from SUBROUTINE PSAP. This 


information does not appear on any Calcomp plots. 
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(NOTE: ONE DATA DECK FROM 
SAP IY CONTAINS ALL 


LOAD 


UCL CARDS 


Figure 4 - 


NAMELIST OPTION 










NAMELIST PICT 
CLORDCASE NO. 2, KODE=) 


CASE 1D CARD 

CLOADCAS: NO. 2) 
NAHELiST PICT 
(LOHDCASE NO. 1/KODE=Z) 








CASES) 
















DISPLACEMENT DATA DECK 
(FROM SAP IV) - 


CASE iD CARD 


GEOMETRY DATA DECK (SRP IY) 


SAMPLE SUBROUTINE PSAP EXECUTION DECK WITH 
IDCASE=1 
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This deck contains the displacement values of the 
nodal points. The deck has one of the following forms, 
depending on the value of KDATA specified in Namelist 
OPTION. 

KDATA = 9 
Calls subroutine DATA9, which reads SAP IV displacement 
data. 


(mete song le ecard, format (15,F10.0), containing 
the number of total load cases (NTLC) that are in the SAP IV 
output displacement data deck, and a scale factor (SCALEF) 
that is used in scaling the displacement data, must be input 
before the displacement deck. This card controls the input 
of SAP IV displacement data through subroutine DATA9. For 
example, a displacement deck from an execution of SAP IV 
that contains the maximum of five load cases would have NTLC 
= 5. The scale factor could be any value desired by the 
user, with default equal to 1. The DATA9 control card for 
the above case would be a5 in card column five and the 


desired scale factor (SCALEF) in card columns 6-15. 


(b) The deck of displacement data obtained from SAP 
IV follows the DATA 9 control card. Since the maxinun 
humber of load cases that are punched out by SAP IV is five, 
the maximum NTLC is five and the actual number aust be 
specified. Note that if NUDISP, NVDISP, and NWDISP are all 
specified as zeroes through default or in Namelist OPTION, 
the displacement data deck and the DATA 9 control card are 
not required. This feature enables preprocessing only of a 


given finite element model by PSAP. 
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(c) The parameter DMAGS in Namelist PICT also 
provides for magnification of displacements. (See section 
III.B.9) 

KDATA = 1 

Calis subroutine DarA1l, which is prepared by the user 
to read displacement data from a program other than SAP IV 
to be plotted. 
KDATA = 3 

Calls subroutine DATA 5, waoich is prepared by the user 
to read displacement data from a program other than SAP IV 
to be plotted. 


Use of KDATA =1 or 5 requires that SUBROUTINE PSAP 
be modified to fit the format that the user's input 
displacenent data specifically requires. A method for 
accomplishing this change will be discussed in parapgraph C. 


of this section. 


The format for Namelist PICT is the same as that 
required for Namelist OPTION. It requires a single card, 
6PICT, followed by the specified parameter cards and the 
mDecabds,, al) Cards beginning in card column 2. This 
NMamelist contains the values needed to specify the type of 
plot that is desired and what information is to be included 
on the flots. A detailed summary of Namelist PICT is 
contained in Reference 2 and is given here for user 
convenience. 

FORTRAN name ~- Default value - Description 
KHORZL - 1 Integer designating the horizontal axes of 
the viewing plane where 1 Xue a= ly and 3=Z,. 


(see Figure 5.) 
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VERT = 


cue — 0. 


foela - 0.9 


foe — 0. 


NEWFR - 


BSCALE - 


ELOISZ = 


XORGN - 


YORGN - 


PSCALE - 


NOTAT - 


med = 0.15 


2 


0 


0 


1 


ORO 


0.0 


0 


KDISP - 0 


Integer designating the vertical axes of 
the viewing plane where 1=X,, 2=Y,,and 3=2Z, . 
(see Figure 5.) 
Angular rotation of model about its X-axis in 
degrees (must be performed third, see Figure 
53) 
Angular rotation of model about its Y-axis in 
degrees (must be performed second, see Figure 
3+) 
Angular rotation of model about its Z-axis in 
degrees (must be performed first, see Figure 
ie) 
1 Frame change before plotting 

(a frame change resets the Y-origin past 
previous plot by XSPACE given in Namelist 
OPTION and resets the X-origin at 0.0) 
0 no frame change before plotting 
1 automatic scaling of plot and computation 
of proper origin location 
2 user-specified origin and scaling 
Maximum dimension desired on completed plot, 
in inches (used for scaling if ISCALE=1) 
X-location of plot origin 

(used if ISCALE=2) 
Y=-locdtion of plot origin 

(used if ISCALE=2) 

Model size reduction factor 

(i.e., PSCALE is equal to actual nodel size 
divided by desired plot size, used 

if ISCALE=2) 

0 no numbering on plots 

1 numbering of grid points 

2 numbering of elements 

Height of integers specified by NOTAT, in 
inches (must be 2 0.07) 


0 plot of undeformed structure 
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plot of deformed structure 


1 

2 exploded plot 

3 displacements represented by vectors 
1 


IDMAG - 2 direct magnification of displacement data 
by DMAGS 
2 scaling of displacement data to a 
maximum value ofDMAGS 

DeAGS — 1.0 Magnification of displacements 
(QESKDISE= 1 or 3) 
Reduction factor of elements 
(L£ KDIS P=2) 

hor axy — 0 1 symmetry about X-Y¥ plane 

KSYMXZ - 0 1 symmetry about X-Z plane 

horMyzZ - 0Q 1 symmetry about Y-Z plane 


Symmetries are perrormed consecutively (i.e., a plate 
quadrant with KSY¥YMXZ and KSYMYZ equal to one would yield a 
complete flate ). 


meex, XXMIN — 1.05420 Locate cutting planes parallel to 
principal planes 

maxX, YYMIN ~- 1.0E~20 Cit 2, i-2)) to fimst plot 

ZZMAX,ZZ2MIN 


momar 9999999 Maximum grid point identification 
humber to be included in the plot 
NDMIN - 0 Minimum grid point identification 
number to be included in the plot 
HEEaak ~ 9999999 Maximum element identification 
number to be included in the plot 
NELMIN - 0 Minimum element identification 
humber to be included in the plot 
KODE - 0 Specifies control option after a 
plot is complete 
0 last plot, exit from program 
1 read another Namelist PICT 
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2 read a new set of displacement 
data (see NOTE 1 below) 
3 read a complete new set of input 
data starting with a title card. 
(NOTE 1: For SAP IV displacement data, KODE=2 signifies 
that the next load case displacements will be assigned to 


the model nodal points.) 


The previous sections describe a complete basic set 
of input data, if KODE =0 in Namelist PICT. For KODE =1, 2, 
or 3, additional sections of the deck must be repeated. The 
deck must end with a Namelist PICT having a value of KODE = 
Jaen | 6 t. An example input data deck and output flots for 
the simple truss problem of Appendix A is found in Appendix 
B. 


Cc. METHOD FOR ALTERING SUBROUTINE PSAP 


In the event a user has geometry and displacement data 
decks from a program other than SAP IV, it is possible to 
plot those decks with PSAP. The subroutine will handle 
rod-like elements, triangular elements, or quadrilateral 
elements when they are input in acceptable format. By 
studying subroutines GEON9 and DATA9, and SUBROUTINE PSAP 
(Appendix D), the necessary sequence of input can be 
determined. PSAP is presently constructed so that a user 
may supply his own routines through the use of subroutines 
GEOM1 or GEOM2 for geometry data and DATA1 or DATAS for 
displacement data. In order to add a subroutine to PSAP, it 
need only be placed after the main calling program in the 
sequence of control cards as discussed in paragraph I1.8.2. 
The essential features for the input of the geometry data 
are the nodal points, with their X,Y,%2 coordinates, and the 


connectivity sequence for the finite element model. The 
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necessary part of displacement data input is the node point 
number, with the U,V,W displacements. Adding a user-prepared 
subroutine through GEOM1 or 2 and DATA? or 5 should prove to 
be a relatively straight-forward task for the user who 
desires to do so. The listing of GEOM9 in Appendix D can be 


used as a guide. 


D. SIGNIFICANT ASPECTS OF SUBROUTINE PSAP 


Often the absence or presence of elements ina 
finite element model cannot be determined from a 
conventional obligue orthographic projection. For example, a 
line element that is coincident with an edge of a triangular 
or quadrilateral element could not be detected. To show 
Clearly each element, PSAP contains an algoritha £Or 
generating exploded oblique orthographic projections. This 
can be a valuable tool in checking the toplogy of an 
analytical model. (i.e., KDISP = 2) 


Zee Portions of Hodels 





The ability to isolate a portion of a model for detailed 
examination is a very useful and desirable asset of a 
preprocessor. SUBROUTINE PSAP has the capability of 
specifying cutting planes (i.e.,XXMIN ,XXMAX ,YYMIN,YYMAX, 
etc.) or maximum and minumum element numbers (i.e., NELMIN, 


NELMAX). Examples of this are shown in Appendix D. 
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The specification of view of a model is done through 
the use of the parameters KVERT,KHORZ, PSI, THETA, and PHI 
fiat are found in Namelist PICT . (Section III.8.9) The 
specific details of how this is accomplished can be found in 
Reference 2. There are a great number on possible 
combinations of the above parameters, and Appendix Cc 
illustrates various combinations that were used on a section 
of a finite element wing model. The Calcomp plots, along 


With the parameters as specified, are found in each figure. 


(a) The most probable source of error could be 
incorrect data deck preparation or deck sequencing. It is 
important to use Figure 2 as a guide while preparing the 
moput for PSAP. 


(6) Other errors nay occur if the arrays 222 and 
DISPD are not dimensioned correctly as discussed in 
Pabagtaph III.8.3. During the execution of the progran 
several manipulations are performed with the two arrays, and 
it is possible for addresses of the data to be lost within 
the IBM 360/67. The 360 System error messages may not 
directly indicate that the dimensions of the arrays have 


been exceeded. 


(c) A most important point to renember in generating 
a sequence of plots is that once a parameter has been 
assigned a value in a namelist, it retains that value until 
it is reassigned. For example, if PLOTSZ is assigned a value 


Sieowomecor the first in a series of plots and it is not 
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redefined in any subsequent Namelists PICT , the value of 
PLOTSZ will retain the value (8.0), as originally specified. 


(dq) It is possible to make any number of errors, 
however, and all of them cannot be anticipated. The error 
messages from the IBM 360/67 System will, in most cases, be 
straight-forward and facilitate easy trouble-shooting. 
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General purpose finite element structural analysis 
programs like SAP IV are significant analytical tools that 
can be used most effectively when coupled with a flexible 
pre~ and postprocessor. The availability of a powerful 
structural analysis coal W2oth Partial pre- and 
postprocessing capability has now been provided for Naval 
Postgraduate School Students. However, at present this 
capability is limited to pre- and postprocessing the data of 
only five of the SAP IV elements. It does not apvear 
feasible, at this time, to add any more capability as far as 
the number of SAP IV elements goes. However, other 
additional improvements to postprocessing at NPS should be 
made. For example, Appendix cC of Reference 2 lists a 
program developed to produce contour plots of stress data 
from finite element models. The program, in all probability, 
could be adapted to SAP IV and the NPS IBM 360/67. This 
represents a possible avenue to further the present 
structural analysis capability for NPS students. SAP EV 
itself can be expanded. There is a great deal of work being 
done in the area of composite materials, and SAP IV 
possesses only very limited orthotropic material capability. 
However, there are subroutines within SAP IV that could be 
modified to reflect more current methods in handling 


ccmposite materials. 


In conclusion, the finite element method in structural 
engineering is the most powerful analysis tool available 
today and needs to be exploited by NPS students in both 
classroom and research work. Programs like SAP IV are 


available at many Navy laboratories, development centers, 
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and aerospace companies throughout the country. The current 
trends in structural analysis are toward computer-aided 
technigues, and exposure to a general purpose package for 
students enrolled in AE 4101 and AE 4102 is a significant 


asset. 
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APPENDIX A 


EXAMPLE PROBLEM USING SAP IV 


The schematic given below illustrates the example 
problem ror which the data on the following page were 


prepared. 


PMS JOM TED 
RSS 
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lL = 100 inches 
B= 10212+06 psi 
ALPHA = 12. 6E-06 psi 


A = 0.5 sq.in. 
Ao =O Sq.in. 
A=) 0575 Sq.in. 
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APPENDIX B 


EXAMPLE PROBLEM USING SUBROUTINE PSAP 


The pin-jointed truss problem in Appendix A is 
as an example in illustrating SUBROUTINE PSAP 
entire computer card deck used to generate Figures 
8 is listed on the following two pages of this 


The printed computer output for this problema is 


also used 

usé. The 
o, 7, and 
appendix. 


shown in 


Figure 9. Note that the displacements of node number 4 have 


been magnified by 1000. 
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Figure = DEFORMED TRUSS MODELS NODEY DISPLACED) 
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Figure a DISPLACEHENTS©NODEY 2? FOR TRUSS MODEL PLOTTED 
AS A VECTOR 
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APPENDIX C 


Searle SUBROUTINE PSAP OUTPUTS FOR A FINITE ELEMENT MODEL OF 
A WING 


The figures contained in this appendix are of a finite 
element wing model developed in. Flight Vehicle Structural 
Amalysis Ii {AE 4102, QTR.IV, 1976). The original wing was 
designed in AE 4274, which is the latter half of the 
subsonic structural design sequence in the Department of 
Aeronautics. The finite element model of the wing was 
developed in AE 4102 as a class project. (The structural 
analysis and design sequences in the the Department of 
Aeronautics are constructed so that they interelate the two 
different areas.) Each student was assigned responsibility 
for a portion of the wing and the student whose design was 


used was the analysis leader. 


Figure 10 shows the axis orientation used in the 
original definition of the wing as designed in Az 4274, 
Figure 11 and 12 show the entire model with over four 
hundred membrane elements and with different views of the 
model specified. In Figures 13 and 14 a single portion of 
the wing is shown with different rotation angles. The 
significant Namelist PICT vaiues used to orient the model 


are included in the figures. 
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WING TIP 





HWING ROOT 


Figure 10 - WING MODEL AXIS ORIENTATION 
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KHORZ=1, KVERT=3 , PHI=70, 
THETA=- 25, PS I=-10 


Fagure 11 °= ORTHOGRAPHIC PROJECTION OF WING MODEL -1 
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=2,KVERT=3, PHI=60, 
THETA=-25, PST 
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MODEL 


ORTHOGRAPHIC PROJECTION OF WING 


Pigure 


56 








KHORZ=2 , KVERT=3, 
PHI=30, THETA=15, 
P§ 1=20, NELMIN=310, 
NELMAX=359 





Fagure 12 = WING MODEL SECTION - ORIENTATION ~-1 
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KHORZ=2 , KVERT=3, 
PHI=60, THETA=-25, PSI=-10, 
NELMIN=310, NELMAX=359 







Figure 14 - HENG MODEL SECTION — ORIENTATION -2 
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APPENDIX D 


SUBROUTINE PSAP LISTING 
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